gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十一章/eg11_3.m

    k=0:0.01:100;
for index=1:10000
p=[2 15 27 k(index)+12 k(index)+1];
r=roots(p);
if max(real(r))>0
break;
end
end
sprintf('系统临界稳定时K值为:K=%7.4f\n',k(index))